Methods for the Study of Particle Production Fluctuations 
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We discuss various measures of net charge (conserved quantities) fluctuations proposed for the 
identification of critical phenomena in heavy ion collisions. We show the dynamical component 
of fluctuations of the net charge can be expressed simply in terms of integrals of two- and single- 
particle densities. We discuss the dependence of the fluctuation observables on detector acceptance, 
detection efficiency and colliding system size and collision centrality. Finally, we present a toy model 
of particle production including charge conservation and resonance production to gauge the effects 
of such resonances and finite acceptance on the net charge fluctuations. 
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I. INTRODUCTION 

The numbers of particles produced in relativistic nu- 
clear collisions differ dramatically from collision to col- 
lision due to the variation of impact parameter, energy 
deposition, baryon stopping and other dynamical effects 
|l], ^, |j| . Such fluctuations can also be influenced by novel 
phenomena such as disoriented chiral condensate pL || or 
the appearance of multiple event classes 0. Even glob- 
ally conserved quantities such as net charge, baryon num- 
ber and strangeness can fluctuate when measured, e.g., 
in a limited rapidity interval. The rapid hadronization of 
a quark-gluon plasma (QGP) can reduce net charge fluc- 
tuations compared to hadronic expectations |) , while 
phase separation can increase net-baryon fluctuations || . 
Fluctuations of conserved quantities are possibly the best 
probes of such novel dynamics, because conservation laws 
limit the degree to which final-state scattering can dissi- 
pate them. 

Many statistical measures have been suggested for 
analyzing particle number fluctuations in experiments 
P, M, hq yj. Although these measures superficially ap- 
pear to be different in nature and unrelated, closer ex- 
amination reveals they are in fact connected. On the 
other hand, each measure exhibits different dependence 
on collision centrality, detector acceptance (rapidity and 
Pt region used to calculate the observable), particle de- 
tection efficiency, and susceptibility to experimental bi- 
ases. The utility of each measure depends on the particle 
species measured and the physical phenomena one wishes 
to extract. For example, "robust" efficiency-independent 
measures are best for observing the correlations between 
neutral and charged particles produced by disoriented 



chiral condensate WIS 



Vc 



Experimental efforts to measure event-by-event fluctu- 
ations have followed two approaches. Many advocate a 
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statistical approach in which fluctuations of particle num- 
bers are characterized by variances, covariances or other 
moments g g, |o[ |ll], [U[ p|, [L6|. These moments can be 
compared to expectations based on thermal-equilibrium 
or other statistical models; any difference can be at- 
tributed to novel dynamics. Others emphasize the im- 
portance of the momentum-dependent correlation func- 
tions, such as the balance function [Til. The correlation- 
function approach has yielded great success in the case 
of identical pion Hanbury-Brown and Twiss correlations. 

In this paper, we discuss relations between correla- 
tion functions and moment measures of net-charge fluc- 
tuations to study the dependence of these measures on 
collision centrality, experimental efficiency and accep- 
tance. We focus initially on the variance Vdyn suggested 
in |Tc| , which is derived from integrals of the single- and 
two-particle distribution functions. Next, we compare 
these measures to alternatives suggested in [0, |n|. Our 
correlation-function based analysis complements a st udy 
by Mrowczynski using a statistical point of view pfj| . 
Specifically, we begin in Sect. II by defining the fluctu- 
ation measure Vd yn in relation to the microscopic corre- 
lation functions. In the next section, we determine the 
scaling properties of Vdyn with system size and, equiva- 
lently, collision centrality. We then introduce alternative 
fluctuation measures, and discuss their relationship with 
Vd yn m Sect. 



[V . A relation between the net-charge fluc- 



tuations, and the balance fun ction intro duced by Bass et 
p7| is presented in Sect. 



al. [y_7j is presented in Sect. [VII]. Sect. |VII]| is devoted 
to a discussion on the robustness of fluctuation observ- 
ables, i.e., whether and how fluctuation measures intro- 
duced in Sect. IV depend on detection efficiency. Finally, 
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we consider and compare, in Sect. IX, the various fluctua- 
tion measures in the context of simple particle production 
models. 



II. NET-CHARGE FLUCTUATIONS AS A 
MEASURE OF 2-PARTICLE CORRELATIONS 



In this section we show that multiplicity fluctuations 
are driven by intrinsic 2-particle correlations. Statistical 
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quantities that we discuss are constructed from the one- 
body and two-body densities, 



dN 

pivn) = -T-, 

CLT] 



d 2 N 
dr]idr] 2 



(1) 



For simplicity, we focus on pseudo-rapidity dependence, 
although results can be generalized to address transverse- 
momentum and azimuthal-angular dependence. Our ap- 
proach and notation in this section follows Ref. (TJ], [2(J . 

To extract statistical information from these micro- 
scopic densities, we use (|l|) to write the multiplicity in 
the rapidity range A77 as 



(AO = / Pl (rj)drj. 



(2) 



Here (N) represents an average of the observable N over 
an event ensemble. Fluctuations of the particle number 
in this rapidity range are determined by integrating the 
two-particle density, 



(N(N-l)) 



Pi(r)\-,r)i)df) X dri2 



(3) 



Ai/ 



The appears on the left side because the integral over 
Pi(j]\-,f]'i) counts the average number of particle pairs in 
the rapidity interval. Note that (N) and (N(N — 1)) 
are the first and second order factorial moments of the 
multiplicity distribution. 

A familiar statistical measure of particle number fluc- 
tuations is the variance, 



V = ((N (N)f). 



(4) 



We can obtain the variance from (|) and (§), since V — 
(N(N - 1)) - (AT) ((AT) - 1). In the absence of particle- 
particle correlations, the two-body density factorizes into 
a product of two one-body densities. In that case, we find 



Ari 



Pi(Vi)Pi(V2)d m d m = (N) 2 . (5) 



The variance is then V = (N), as expected since the 
number of particles produced in a sequence of indepen- 
dent events follows Poisson statistics plj . Note that 
the relative uncertainty in the mean number (AT) is 
VV I (N) = 1/ yj (N) for this case. Observe that the par- 
ticle number in a grand canonical ensemble in thermal 
equilibrium follows Poisson statistics. 

Information on net-charge fluctuations is contained in 
the two-body density for distinct particles with opposite 
charges. We determine these fluctuations from 



(N a Nfs) 



P2(Va,rii3)dr) a dr)/3, 



(6) 



where a and [3 label the particle species. In a statisti- 
cal framework, this average is related to the two-particle 
covariance, 



V a p = (NaNp) - (N a )(Np). 



(7) 



The covariance vanishes if there are no correlations 
between the species a and (3, since p2(?7cn?7/3) = 

Pi(Va)pi{vp)- 
Following pi 



20 we define the robust variance, 



Rcta 



(AT) 2 ' 



and the robust covariance, 



R, 



■a(3 



(8) 



(9) 



for particle species a and (3. These quantities have 
the same sensitivity to fluctuations as the variance (^) 
and covariance (0) but have three significant advantages. 
First, these quantities vanish for V = (N) and V a $ = 
0, so that they measure the deviation from Poisson- 
statistical behavior. Second - and of greater practical 
importance - the ratios (||) and (||) are 'robust' in that 
they are independent of experimental efficiency, To see 
why Eq. (^J) is robust, let the probability of detecting each 
charged particle be e and the probability of missing it be 
1 — e. For a binomial distribution the average number of 
measured particles is (Af) oxp = e(N) while the average 
square is (N 

2\ 

exp 



K X p = (iV 2 ) cxp ~(iV) c 2 xp = e 
so that V c 



,., p = e 2 (N 2 ) + e(l - e)(N). The variance 
2 c ^ = e 2 ((N 2 )-(N) 2 ) + e(l-e)(N), 



cxp 



(AO exp = e 2 {V- (AT)). We then find 



noxp n 



(10) 



independent of e; the proof that (|9|) is robust is similar. 
The ratios (||) and (j^) are strictly robust only if the ef- 
ficiency e is independent of multiplicity. We discuss this 
point in more detail in section VIII. 

Third, R a p are directly related to the particle corre- 
lations. For a ^ /3, we combine (|^), (^|), (0) and ([)]) to 
obtain 



R a p = 



J Av P2(Va,Vi3)dri a drip 
J Av Pi(v a )dv a f An pi(r)p)dr)0 



(11) 



one can check that ( |11[ ) also holds for a = ft. As in an 
HBT analysis, we define a correlation function C by 

P2(vuV2) = pi(m)pi(m)[i + c(m,m)}, (12) 

so that ( pi] ) yields 

J Av Pl{V a )Pl(Vfi)C a t3{Va,Vt3)dT] a dl] l 3 



Rafi — 



(13) 



We use this result to illustrate how to extract microscopic 
information on the rapidity range of correlations from the 
A77 dependence of R a j3 in sec. VII. 



3 



To study net-charge fluctuations, one can measure the 

robust covariance for charged hadrons R^ On the other 

hand, it would be better to isolate the potentially inter- 
esting net-charge fluctuations from factors that cause the 
numbers of positive and negative hadrons to fluctuate to- 
gether, such as variations in energy deposition or collision 
volume. Toward that end, we consider dynamic charge 
observable defined as the linear combination 



Vdyn — 



(14) 



Ratio fluctuations considered by Jeon and Koch 
an alternative 



see sec. IV 



This combination vanishes 
when the positive and negative hadrons fluctuate simul- 
taneously, since all the R a p are then the same. We also 



see that Vd yn is both robust (see sec. |VIIlD and straight- 
forwardly related to the microscopic correlators (|l|), as 
are the R a p- We find an alternative expression for Vd yn 
in terms of 



( N + 




(15) 



where N+ and N— are respectively the multiplicities of 
positive and negative hadrons. In the limit of indepen- 
dent particle production, v becomes 



1 



1 



v stat 



(N+) (JV_) • 
The dynamic charge observable is the difference, 



Vdyn 



V - V s tat, 



(16) 



(17) 



as we see by expanding the square in (|15|). Observe that 
Vdyn is nonzero when net-charge fluctuations are corre- 
lated (non-Poissonian) . Furthermore, eqs. ( |l5[jl7| ) are 
more useful than ( |l4| ) for extracting correlations from 
numerical data since the net-charge fluctuations are typ- 
ically smaller than the fluctuations of the total number 
of hadrons. 

We examine the scaling properties of the Vd yn variance 
with collision system size in the next section. 



III. 



SCALING OF 



l',lu 



WITH SYSTEM SIZE 



AND COLLISION CENTRALITY IN A+A 
COLLISIONS 

We now study the scaling of the observables C a p, R a p, 
and Vdym with collision centrality, target and projectile 
mass. For concreteness, we assume that nuclear colli- 
sions are a superposition of independent nucleon-nucleon 
(NN) sub-collisions and neglect the rescattering of the 
hadrons. These assumptions imply that charged-particle 
pairs can be correlated only if produced in the same sub- 
collision. We expect the contribution to the two-body 
density from these related pairs to grow linearly with the 



number of sub-collisions M . Related pairs will be diluted 
by random pairs. The A A densities are 



92 



pi aa (v) 

AA (Vi,V2) 



Mpi NN ( V ), 
Mp 2 NN ( m , m ) 

NN i 



+M(M l) Pl NN ( m ) Pl NN ( m ). 



(18) 
(19) 



The first term of ([L9|) describes the related pairs while 
the second accounts for the M(M — 1) random pairs. 
These expressions apply generally to particle production 
from M sources; we focus on the independent collision 
model for simplicity. We apply these considerations to 
more realistic models at the end of this section. 

To compute the correlation function, we substitute the 
AA densities Eqs. in ([jj|) to find 



c A £(vi,m) 



M 



(20) 



For independent sub-collisions and in the absence of 
rescattering, we therefore expect the A A correlation func- 
tion to have the same rapidity dependence as in pp col- 
lisions, with an overall scale that is reduced by a factor 

Before turning to realistic experiments, we consider 
for the moment a collision with a fixed number of sub- 
collisions. The statistical observables then satisfy 



and 



nAA 

n aj3 



v dyn {AA) 



nNN 
M 



VdyniPP) 
M 



(21) 



(22) 



We see that all quantities scale as M _1 . 

More realistically, suppose that one specifies a cen- 
tralility range by measuring the total charge multiplicity, 
the zero degree energy, or some analogous global observ- 
able. The number of sub-collisions will then fluctuate, 
adding to the variance and covariance of particle num- 
bers and changing (pil"|). Specifically, the fluctuations of 
M contribute a term (N a ) (N p )((M 2 ) - (MS 2 ) to the vari- 
ance and covariance, V n and V a p , so that (0) and M) give 



jdNN 
r>AA _ n af] 
U a(3 — 



(M) 



(M 2 ) - (My 



(23) 



See the appendix for a full derivation. We remark that 
these M fluctuations are essentially equivalent to the 
"volume fluctuations" discussed in a local equilibrium 
framework jj], |§| . 

On the other hand, random changes in the number of 
independent sub-collisions can change the total number 
of particles but not the net charge, so that (53) is effec- 
tively unchanged. We find 



Vdyn(AA) = R ++ + - 2R + . 



VdyniPP) 

(M) ■ 



(24) 
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The contributions from sub-collision or volume fluctu- 
ations are the same for all a and (3, so that ( |f4| ) im- 
plies that this contribution does not affect Vdyn- 

The 

second term in (|2l|) is of order f / (M) and comparable 
to the first, since ISR and FNAL experiments suggest 
that i?£f ~ ~ R±?/2, each of order unity in 

Atj = 1 - 2 at RHIC p), §(J. 

We now extend these considerations to the wounded 
nucleon model, which succssfully describes many global 
features in SPS and AGS experiments. There, one as- 
sumes that only the first sub-collision of each nucleon 
drives particle production and neglects all subsequent 
interactions (22). Since ( |l8| ) and ( |f9| ) formally describe 
particle production from M independent sources, we can 
adapt © and (jl|) to the wounded nucleon scenario 
by replacing the number of sub-collisions M with the 
number of participant nucleons, M. We must also re- 
place the densities p\ NN and p2 NN in ( p"8| ) and ( |f9| ) with 
coefficients p\ and p\ that describe the production per 
participant. Observe that nucleons are counted as par- 
ticpants if they interact at least once and that there are 
two participants per XX collision. 

Results of the form (^3|) and (|24|) then follow from 
the wounded nucleon model if we replace M with one- 
half the number of participants M. . The average number 
of participants at impact parameter b for a symmetric 
AA collision is (M(b)) = 2 J dsT{s){l - e^™ 3 ^-')}, 
where T(b) — J p(z, b)dz is the familiar nuclear thick- 
ness function and p is the nuclear density. By com- 
parison, the number of sub-collisions is (M(b)) — 
cjnn J dsT(s)T(b — s). We remark that both wounded- 
nucleon and independent-collision approximations imply 
that the total multiplicity of pions N n scales as the re- 
spective number (participants or subcollisions) . There- 
fore, both models imply Vdyn °c X" 1 , albeit with different 
coefficients. 

We point out that particle production at RHIC energy 
has contributions from soft interactions, which scale as 
the number of participants, and hard processes, which 
scale as the number of subcollisions |^3|, |24| . In this case 
the scaling of Vdyn with N„ can be more complex. Fur- 
thermore, final state scattering effects can certainly mod- 
ify this scaling. 



IV. ALTERNATIVE MEASURES OF 
FLUCTUATIONS 

In this section we consider the connection between the 
variance Vdyn and other fluctuations measures. We dis- 
cuss some of the merits and problems associated with 
each observable. 



tical considerations. It consists of the difference be- 
tween the mean of particle production variances calcu- 
lated event-by-event and the variance calculated over the 
entire dataset. Consider x an observable of interest e.g. 
the net charge of produced particles. The inclusive mean 
of x (i.e. average over all particles an events) is noted x. 
Deviation from the inclusive mean are noted Air = x — x. 
By construction, one has A x = 0. Th e root mean square 
(RMS) deviation is Ax 2 = (x — x) 2 . To investigate the 
dynamics, one determines how the event-wise net value 
of u x", defined as X — J2i x i' changes event by event. 
One defines AX = X — Nx as the event deviation from 
the inclusive mean (with N being the number of particles 
in the given event). By construction, its event average 
(AX) vanishes, whereas (AX 2 ) does not. The $ measure 
is defined as |l(| 



(25) 



For a system with particles of charge q+ and g_, the 
inclusive standard deviation is: 



Ax 2 = (q+ - q_)- 



(N) 2 



(26) 



The magnitude of (AX 2 ) is determined by both statis- 
tical and dynamic fluctuations. Defining Q as the net 
charge of an event, one has AX = Q — N(Q) / (N), from 
which one finds indeed (AX) = 0. The average (AX 2 ) 
is however non-zero. One finds 



(AX 2 ) = (q + -q_ 



i2 (X + ) 2 (X_) 2 



(NY 



(N 2 + ) (X + ) 2 

(X+) 2 
(Nl) - (X_) 2 



-2 



so that 

$ = (q+-q 



(AM 2 
(N + N_) - (N + )(N_) 

(N+)(N-) 
(X+)(X_) /(X 2 )-(X+) 2 



(27) 



-2 



(JV)a/a V (jv+) 2 
(X 2 ) - (X_) 2 

(X + X_)-(X + )(X_)\ 1/2 



(N+)(N-) 



(N+HN-) 
(N) 2 



1/2* 



(28) 



A. ^-Measure 

The $ measure of the net charge fluctuation was in- 
troduced by Mrowczynski and is based on statis- 



Examination of eqs. (p6|), (p7|), and ( psj ) reveals that they 
can, in fact, be expressed as the v and v sta t variances as 
follows (as also reported by Mrowczynski pCj): 



(AX 2 



(AO 3 



(29) 
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= ( g+ -g-) 2 V ( y )3 



^stai (30) 



so one can express $ as 

2(N+)(N_) ( flT [U^\ 



$ = 



(N) 



(N) V W J 



(31) 



In general, the dynamic component of the fluc- 
tuations is much smaller than the statistical com- 
p oncnt: v dyn « v sta t implying y/y]{N) - 

y/vstat/ (N) = y/u stat /(N)(y/l + Vdyn/ ^stat - 1) ~ 

Vdvn (2 \/ Vstat (N) ) ~ 1 ■ Substituting the value of v sta t 
given by Eq.fLq, the above expression can thus be ap- 
proximated by: 



(JV+) 3 / 2 (A_) 3 / 2 

$ ~ — w 2 — Vdyn ' 



(32) 



One thus finds that indeed the measure is determined 
(mostly) by the dynamical fluctuations of the system, i.e. 
by the particle correlations implicit in the sum R++ + 

i? 2i? + . 

Eq. (|3^) further simplifies, as follows, for cases where 
(N + ) = (N_): 



(33) 



Given, as we discussed in Sect. [II, that the variance v dyn 
should vary inversely to the multiplicity of charge parti- 
cles in the limit of independent particle collisions and 
absence of rescattering of the secondaries, one should ex- 
pect that <!> w Vd yn ,pp/& in that limit, and independent 
of the collision centrality if the collision dynamic do not 
vary with collision centrality. Note however one must 
exercise caution while comparing $ measured by exper- 
iments with different acceptances (See Sect. [Vj for de- 
tails). Note finally that unlike Vdym the $ measure is 
a non robust observable given it explicitly depends on 
the detection efficiency of positive and negative particles 
through the factor (N + ) and (iV_) as we shall discuss in 



more detail in section VIII. 



B. Particle Ratios 

Another approach advocated in ref. || focuses on the 
variance of the ratio of positive and negative particle mul- 
tiplicities, R = (N + ) /(N-). As shown in ref. the fluc- 
tuations of the ratio offer the advantage that "volume" 
fluctuation effects cancel to first order. This is also true 



for Vdyn (see sec. Ill) and $ |10 



For small fluctuations, the variance of the ratio can be 
related to the charge variance v ([[5]) . A small fluctuation 

of R= (N + )/{N-) satisfies 



AR 
~R~ 



A7V + 



A7V_ 



(34) 



so that 
(AR 2 ) 



{AN 2 



(AN 2 :) _ 2 (AN+AN-) ^ (35) 



(R) 2 (N + ) 2 (N.) 2 (N+)(N-) 

Expanding the square in (|l5|), we see that 

{AR 2 ) = (R) 2 v (36) 

Observe that neither v nor {AR 2 ) are robust. Also, note 
that this equivalence holds only when (AN^) 1 / 2 « 
(N±); an approximation, which breaks down at small 
multiplicities. Problems with these quantities for small 
multiplicities are discussed in |l5fl . The D measure used 
by Jeon and Koch B, 

D = (N) (AR 2 ) = (N+ + N_){R) 2 v, (37) 

is also efficiency dependent. 



C. Reduced Variance 

Lastly, we consider the reduced variance ujq used by 
authors f| §, 0, [D| §||. If we write N = N + + N_ and 
Q = N + — N- , then the reduced variance is 



(AQ 2 

(N) 



As before, we expand the square to find 



(AN 2 ) + (AN 2 ) - 2(AN+AN. 



(N+) + (JV_) 



(38) 



(39) 



This ratio is unity for Poissonian fluctuations or for a 
thermal ensemble in chemical equilibrium; any measured 
multiplicity dependence would be interesting. In terms 
of robust ratios, we obtain 



ujq = 1 



(N, 



(N) 



-R 



++ 



(N-l 
(N) 



R. 



, (N + )(N. 
' (N) 



(40) 

Generally, this quantity has a complicated dependence 
on the correlators R a p. However, for (N + ) w (iV_), the 
above expression reduces to 



f 



<A 4 



N-) 



(41) 



indicating that this quantity has the same efficiency de- 
pendence as the total number of charged particles. 

We note, in closing this section, that the reduced vari- 
ance, u>q, unlike Vdyn, and <&, has an explicit dependence 
on collision volume fluctuations, as given by the following 
expression. 



UQ — + 



((N + ) (AM) 2 <AV 2 } 



(42) 



(N + ) + (N.) (V) 2 

where ujQy corresponds to the reduced variance at fixed 
volume, while (V), and (AV) 2 are respectively the mean 
and variance of the collision volume. The importance of 
volume fluctuations was pointed out by Jeon and Koch 
p6[ . Following their work, it is straightforward to show 
that Vdyn, and $ are independent of volume fluctuations. 
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V. CHARGE CONSERVATION EFFECTS 

The total charge of the system is fixed due to the 
charge conservation. It implies some "trivial" correla- 
tion in particle production regardless of other dynamical 
effects. As such, it only affects the two-particle density 

p H (?7 + ,7/_). We proceed to study the effect of charge 

conservation on the net charge fluctuation by calculat- 
ing the correlation function C_| (r)+,r]-) as a function 

of single and two-particle density expressed in terms of 
probability distributions of positive and negative parti- 
cle in order to emphasize the role of charge conservation. 
One writes, for fixed number N± of positive and negative 
particles: 

P±(v±) = N ± P±(r,±) (43) 
p + _(77+,r;_) = N.P + .(r] + ,r)-) (44) 
+{N-N + -N-)P + (r i+ )P-(j l J) 

P±{ r l±) are probabilities to find one + or - particle at 

rapidity rj±. Pj is the probability to find one positive 

particle and one negative particles at rapidities rj + and 
rj- respectively. iV_ and N + are respectively the total 
number of negative and positive particles produced (over 
47T solid angle) by a collision. By virtue of charge con- 
servation, and given the total charge Q > 0, one has 
Q = N+ - iV_ , and N + > 7V_ . The first term of Eq. (@) 
accounts for correlations between positive and negative 
particles. As there are N- +- pairs created, one has a 

contribution iV_P H The second term arises because 

there are N+N- — A_ ways to pair the uncorrelated +- 
particles. In general, at a fixed impact parameter (or 
number of nn collisions), the multiplicities N- and N + 
shall fluctuate event-by-event. One must then average 
over such fluctuations and rewrite the above expression 
as 

P±{V±) = (N±)^P±(v±) (45) 
P+-(V+,V-) = (N-)^P+-(v+,V-) + (46) 
((N_N + ) i7r + (N.) 4 „)P + (r, + )P^) 

where the notation (0)i n represents an average taken 
over 4-7T acceptance. In the absence of dynamical corre- 
lations, and by virtue of charge conservation, one has 

(N 2 _) 4n - (N.) 2 ^ = (iV_) 47r (47) 
(JV_JV+>4 ir = (iV_)L + (N_) 47T - {N_)^Q 

The correlation function C_| (r;+ , r\- ) can then be calcu- 
lated and written as 



1 



p+-(v+,v-) 
p+{v+)p-{v-) 

1 P+_(r ?+ ,7 ? _) 
(N+) 47r P + ( V+ )P.(r,.) 



(48) 



This result is fairly generic and includes the possibility 
of dynamical spatial (or rapidity) correlations between 



the particles of a created pair. Neglecting such a cor- 
relation however, and for the purpose of evaluating the 

role of charge conservation alone, one sets P_| = P + P_. 

One then finds that charge conservation implies: 



1 



(N. 



+ /47T 



(AO 4. 



(49) 



where (N)^ stands for the mean total number of charged 
particles produced in the event. Obviously, at large mul- 
tiplicities one can neglect the difference between N + and 
N/2. 

The correlator P_| is obtained by integration (see 

Eq. (|TT|)) of C-| (?7 + ,77_) over the experimental accep- 
tance. Given that C + _ (ry + , rj_) is actually independent of 

r]±, -R-i is independent of the experimental acceptance. 

One thus finds that the charge conservation contribution 
to Vdyn amounts to 



Avdyn — — 



(50) 



It is independent of the experimental acceptance, and 
only determined by the total charge particle multiplicity 
at a given impact parameter. 

We emphasize that Vdyn 7^ for a 4-7T acceptance be- 
cause charge conservation imposes a correlation on the 

system. The total given by ([l5]) is strictly zero 

when all particles are detected. However, (|l^) implies 
that Vstat 7^ in this case, since the Poisson distri- 
butions used to calculate v sta t do not incorporate a 
global charge conservation constraint. It follows that 
fdvn = v — v s tat — > —v s tat for a An acceptance, as seen in 
(pp|). This estimate of the effect of charge conservation 
is in agreement with a correction reported in p7| . Note 
however that the correction is additive not multiplicative 
as stated in [p7l. 



VI. RAPIDITY DEPENDENCE OF 
FLUCTUATIONS AND DETECTOR 
ACCEPTANCE 

Measuring the dependence of R a p and Vd yn on the ra- 
pidity window A?y can yield information on the rapidity 
range of correlations as well as their magnitude. Informa- 
tion on the rapidity dependence of R a /3 is also needed to 
compare data from experiments with different geometric 
acceptance. The microscopic correlations themselves can 
and indeed must be determined from balance function 
and similar measurements Ji"7[ ; such experiments have 
different practical issues. We relate Vd yn and balance 
function measurements in the next section. 

To exhibit the rapidity dependence of Rap, we as- 
sume that p\ are rj independent and that C — 
C a fl(0)exp{-(?7i - 77 2 ) 2 /2ct 2 }. ISR and FNAL data 
|l9| , |20f1 show that charged particle correlation are func- 
tions of the relative rapidity r/i — r\2 with only a weak 
dependence on the average rapidity of the pair. Data 
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FIG. 1: Rapidity dependence of the robust covariance R+- 
assuming a Gaussian correlation function of width a. 



can be roughly characterized as Gaussian near midrapid- 
ity. Using (O) we find 



Ra/3 



{V^erf(»- (l-e-* 2 )}, (51) 



where x = \f2Ar\ja. The function R a p is shown as a 
function of A77 in Fig. 1. ISR and FNAL data suggest 
that the rapidity range of correlations is roughly from 1 
to 2 rapidity units. 

Both R and the microscopic correlator C depend on 
the value C a p(fS) at rji = r\ 2 and the rapidity range of 
correlations, a. Equation ( |5l| ) carries the same value - 
and caveats - as does the Gaussian parameterization of 
HBT correlations. The range a depends on the dynam- 
ics and may vary with centrality as well as target and 
projectile mass. 

One must account for this rapidity dependence when 
comparing experiments of different geometrical accep- 
tance. We estimate, for instance, that the difference be- 
tween of the fluctuations measured by the STAR (\rj\ < 
1.5,A0 = 2vr) and PHENIX {\q\ < O.35,A0 = tt/2) ex- 
periments to be roughly ~ 10% for a ~ 1 — 2. While 
this is a rather small correction, we emphasize that the 
experiments should measure the rapidity dependence. In 
general a can differ from pp to AA collisions and, more- 
over, is expected to depend on centrality. 



VII. RELATION BETWEEN THE BALANCE 
FUNCTION AND v dyn 

The balance function was proposed by Bass et al. JTtJ 
as a technique to study the dynamics of hadronization 
in relativistic heavy ion collisions. The idea is that the 
rapidity range of correlations is changed when a colli- 
sions forms quark gluon plasma. Specifically, charged 
hadrons form late in the reaction, after hadronization, 
resulting in shorter-ranged correlations in rapidity space 
for charge/anti-charge pairs than expected in the absence 
of plasma. 

The balance function as defined by Bass et al. jl7| is 
written (here again focusing, without loss of generality 



on the rapidity dependence) 

1 



B(A77 2 |A?7i 



^ {£>(-, A/7 2 |+, At?!) 

-D(+,A»fc|+,A»7i) 
+£>(+, A772K A^) 
-D(-,A m \-,A m )}, 



(52) 



where 



D(b, Ar) 2 \a, A Vl 



pi+Ar;i/2 r) )2 +Ar ?2 /2 
Jt/jl-Atji/2 J))2-Aj )2 /2 



rri2+Ar) 1 /2 , , , 

Jm-A m /2 dVaPliVa) 

(53) 

The ratio D(b, Ar/2\a, Ar]i) is essentially a conditional 
probability for finding a number of particles of type b 
in the phase space bin A772 of centroid 772 given the pres- 
ence of particles of type a in the phase space bin A771 of 
centroid 771, i.e. 

A^(6,A77 2 ;a,Ar ?1 ) 



D(b,Arj 2 \a, A771) = 



(54) 



N(a,A m ) 

The bins need not overlap. Experimentally, evaluations 
of the balance function can be restricted to a determina- 
tion of the correlation of particle a and & as a function 
of their relative rapidity Ary. In this case, particle a can 
be anywhere within the full detector acceptance Y, and 
particle & is at a rapidity A77 relative to a. This leads to a 
one-dimensional balance function, B(Arj\Y), defined as: 

B(Ar)\Y) = -{D(-,Arj\+,Y) 



-D(+, Ar)\+,Y) 
-D(+,Arj\-,Y) 
-D(-,A V \-,Y)}, 



(55) 



To understand this expression better, observe that for a 
sufficiently narrow bin A77 we can write 

An [ Y ' 2 

D(b,Arj\a,Y)^—L d Va p 2 (r, a ,r)), (56) 

\Ma) J-Y/2 

where (N a ) is the number in the full domain —Y/2 < rj < 
Y/2. For a boost invariant system the pair correlation 
function C is a function only of the rapidity difference, 
so that this integral is essentially C averaged over the 
system volume, plus a constant term that cancels in (|5^). 

The integral of this function over the entire acceptance 
Y is noted B(Y\Y). By virtue of Eq. (|54|), it amounts to: 

B(Y\Y) = [ dArjB{Ar]\Y) 
Jo 

(N + N_) Y 



(N+)y 
(N + N^) Y 

(JV_)y 

(N + (N + - l))y 

(N+)y 
(N-(N--1)) y 



(57) 
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The four terms of this equation are part of the expression 
of the correlators R a h given in Eqs. d||l]). The integral 
B(Y\Y) can thus be re- written: 



B(Y\Y) = -{i? + _(A_)+i? + _(A + ) 
-R++{N + ) - R—(N_)} 



(58) 



which establishes a relationships between the integral, 
B(Y\Y), of the Balance function, and the correlators 

R++, R--, and i? H 

At RHIC, one observes that (JV_) » (N+) = (N)/2 
near central rapidities in Au+Au collisions p8[ ]. The 
above expression simplifies 



B{Y\Y) 



- - R — } 



(59) 



The integral, B(Y\Y), of the balance function, B(Ay\Y), 
is thus indeed proportional to the variance Vdyn and the 
total multiplicity (N) when (iV_) » (iV + >. 



VIII. FINITE RECONSTRUCTION 
EFFICIENCY EFFECTS 

We consider the effect of finite reconstruction efficiency 
on measurements of fluctuations studied as a function 
of collision centrality. We assume the centrality is ex- 
perimentally determined based on the total multiplicity 
of charged particles detected in a reference acceptance, 
f2jvf whereas the multiplicity fluctuations of interest are 
measured in a fiducial acceptance f2jv- We account for 
the finite detection efficiency, in a given acceptance, Q a , 
by introducing a detector response function Po(n a \ N a ) 
expressing the probability of detecting a multiplicity n a 
given a produced multiplicity N a . In general, Pd (n a \N a ) 
shall account for finite efficiency effects as well as mea- 
surements of ghost tracks. We shall calculate, quite gen- 
erally, moments, Mfc jQ , and factorial moments, Fk >a , of 
the particle multiplicity distribution defined respectively 



Fk,a 



i 

Aei 

(N a (N a -l)(N a -k)) 
-^^AT a (]V a -i)...(JV Q -fc), 

1 * PA) 



(60) 



IS H a 



where N ev is the number of events studied. The mean 
= Afi, Q and the variance, V = (SN^) = M 2: a — 
Here we will restrict our calculation to these lowest 
moments, but the calculation can easily be generalized to 
higher moments. 

We shall use lower case letter (e.g. mk, a ) to distinguish 
measured moments from the intrinsic or actual moment 



of the produced particles (i.e. that one wishes to infer) 
represented with capital letters (e.g. Mk >a ). 

We assume that moments of the multiplicity distribu- 
tions are measured as a function of the collision centrality 
estimated based on the total multiplicity, m, measured in 
the reference acceptance. The moments can then be ex- 
pressed (neglecting for simplicity the particle type label 

a): 



m k 



\ n k P(n\m), 



(61) 



where the sum is taken over all relevant multiplicities, 
and P(n\m) is the probability to measure "n" given the 
centrality estimator W . We emphasize that both "n" 
and "m" are influenced by the finite efficiency of the de- 
tector. We in fact seek to extract the intrinsic moments 
of the particle production 



M k = N k P{N\M), 



(62) 



N=0 



where P(N\M) is the probability "N" particles are pro- 
duced at a given centrality "M" . The measured distribu- 
tion P{n\m) can be expressed as a function of the intrin- 
sic distribution as follows 

P(n\m) = ^2 PD{n\N)P(N\M)P D {M\m), (63) 

N,M 

with the sum extending over all relevant produced mul- 
tiplicities N and M. The factor Pjj(M\m) corresponds 
to the probability of having a produced multiplicity AI 
given the measured value m. It is evaluated using Bayes 
rule: 



P D {M\m) 



P D {m\M)P{M) 
P(m) ' 



(64) 



where P(M) and P(m) are respectively the probability of 
the produced, M, and measured, m, multiplicities. The 
measured probability distribution is thus 



P(n\m) 



1 



£P(n|JV) 



P(m) 

xP(N\M)P(m\M)P(M). 



(65) 



Measured moments can be calculated as function of 
the intrinsic (produced) moment by inserting the above 
expression in (|62|). Introducing for convenience the func- 
tions h s (N) and g s (M) defined as follows: 

h s (N) = ^VP(n|iV) (66) 

n 

g s (M) = ^2h s (N)P(N\M) (67) 

N 

one finds a general expression for the moments as follows: 
<™fc> = TJT^T E P(m\M)P(M)g k (M), (68) 



P(m) 
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Assuming P(n\N) can be appropriately approximated by 
a binomial distribution, the above expressions can be 
readily simplified. The moments h s (N) yield 



h^N) = eN 

h 2 (N) = e 2 N 2 + e(l - e)N, 



(69) 



where e n is the detection efficiency achieved in the mea- 
surement of "n" . Substituting these quantities in Eq. j67| ) 
leads to 



9i(M) 
g 2 (M) 



e n (N) 
el{N 2 ) 



t (l-e n )(N). 



(70) 



The first and second moments, are thus in general 



P(m) 



M 



(n 2 ) = -L_\2 P D(m\M)P(M) 



M 

-2 / Ar2 



x(e 2 n {N 2 }+e n (l-e n )(N)) 



(71) 



with the moments (-/V s ) evaluated at a fixed value of M. 
Clearly, the measured moments are determined by the 
intrinsic moments smeared over the response function of 
the multiplicity, M . Assuming the efficiency of the total 
multiplicity detection process is near unity, one can ap- 
proximate the response function P£>{m\M) with a delta 
function 8 m .M, and the above expressions simplifies as 
follows 



„2\ 



e n (N) 
el(N 2 ) 



i(l-en)(N). 



(72) 



We show in appendix the above results holds for fi- 
nite efficiency as long as "n" has a linear dependence on 
the total multiplicity "m" over the range of the response 
function P D (m\M). 

We now proceed to use these for the calculation of 



the various fluctuation measures introduced in Sect. [TV. 
We use sub-indices "+", "Q" , and "CH" to denote 
positively and negatively charged particles, net charge, 
and total charge particle multiplicity, respectively. We 
use overlined symbols to re pres ent the intrinsic measures. 
We find using Eqs. F|, ||, gl 



1 — e± + e±u)± 
l-e±+ e±uj Q 
1- e±+ e±uj C H 

3/2 3/2 



(73) 



The above fluctuation measures display an explicit de- 
pendence on the charged particle detection efficiencies e± 
or the total efficiency e. The $ observable, in particular, 
has a non trivial dependence on the detection efficiencies 
of positively and negatively charged particles. This de- 
pendence however simplifies to a single factor, e, if the 



positive, negative, and global efficiencies are equal (i.e. 
e + = e_ = e). By contrast, one finds that the dynamic 
variance Vdyn — v<iyn i-e. it is independent of the de- 
tection efficiencies, and is thus, in that sense, a robust 
observable. Note that this conclusion remains strictly 
correct as long as the Gaussian approximation is valid. 
See the Appendix for a discussion of the Gaussian ap- 
proximation. 



IX. SIMPLE PRODUCTION MODELS 
A. Poissonian Particle Production 

We first consider a multi-particle production model 
where no correlation are involved. Specifically, we as- 
sume that on average, particle species, i, are produced 
in fixed fractions of the total particle production. Wc 
consider cases where the fluctuation measures are evalu- 
ated over kinematic ranges that might be identical (case 
A) or smaller (case B) than the kinematic range used to 
calculate the total (charge) particle production. 

The probability to produce species, i, with multiplici- 
ties, Ni, is evaluated with a multinomial distribution. In 
general, one has 



P(N 1 ,N 2 ,---,N k \M) 



1 



n 



N,, 



J a 



(74) 



1 



In case A, one shall have M = N a and fa 
whereas in case B, M > N a , and f a < 1. 

The multiplicity moments, and variance are calculated 
at fixed total multiplicity, M, assumed to be representa- 
tive of the collision impact parameter: 



(N a ), 

(NaNp), 

(N a (N a -l)) r 



f a M 

f a M + f 2 M(M - 
M(M - l)f a f p 
f 2 M(M - 1) 



1) 



V a = Mf a (l-f a ) 



(75) 



Consider now the specific case of net charge fluctua- 
tions with the index a taking values + and — . One has 



in case "B" 



V Q = 

= 

U c h — 
V+- = 

$ = 



1 



(A 



/- - (/-+ 
-f-) 2 



f-f) 



U + f- 

1 - (/+ - /-) 

u - f- 

Mf+f- 



0. 



(76) 



Case A is easily calculated from the above by setting 
/+ -/- = !■ 
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The coefficients f± can be experimentally determined. 
It is thus straightforward to determine the normalized 
variances expected for particle independent production 
and compare with measured values to seek for the pres- 
ence of sub or super Poissonian fluctuations. Note ad- 
ditionally that both the Vd yn and <E> variables have null 
expectation values irrespective of the fraction of the frac- 
tions f±. They thus constitute a more reliable measure 
of the dynamic fluctuations. 



B. Simple Resonance Production Model 

Two-particle correlations are determined by a host of 
phenomena such as collective (flow) effects, production of 
resonances, jet production, Fermi/Bose statistics, as well 
as intrinsic phenomena related to the underlying collision 
dynamics. Here we examine the role of resonance decays 
(e.g. p°, A°) on measurements of the net charge fluctua- 
tions. We show that the production of neutral resonances 
that decay into pairs of positively and negatively charged 
particles produce an effective dynamical correlation. 

We formulate a simple toy model where we include only 
three types of particles: ir + , ir~ , and p° . The p° shall be 
viewed as a generic neutral resonance, which decays into 
7r + and 7r~ . Obviously, this is an oversimplification of the 
problem and a fuller treatment shall account for other 
species, all relevant resonances, and the finite acceptance 
of the detection apparatus. 

We consider the n + , ir~ , and p° to be produced inde- 
pendently (neglecting Bose effects) at freeze out in rela- 
tive fractions f\, f 2 , and f 3 respectively, and model the 
multiplicity production according to a multinomial dis- 
tribution (as in the previous section) . The probability of 
producing n\ 7r + , n 2 7r~, and n 3 p° is expressed 



P(ni,n 2 ,n 3 ;N) 



Nl 



ni\n 2 \n 3 \ 



Jl J2 J3 ■ 



(77) 



Given our assumption that all p° decay into a pair 7r + and 
7r~, the probability of measuring n + positive n_ negative 
particles respectively can be written 

P(n+,n_;iV) = £ P( m,n 2 ,n 3 ;N) 

n 1 ,n 2 ,n 3 

^ &7l+,7l\ +713*^71— ,712+713 ■ V^J 

One then writes the moment generating function of the 
probability P(n+,n_; N) as 

G(t + ,t_;N) = (p ie *+ +p 2 e t - +p 3 e t + +t -) N , (79) 

which one uses to computes the moments of the pion 
multiplicity distributions. One finds: 

(N+) = N(f 1 + f 3 ) 

(N-) = N(f 2 + f 3 ) 

(N+(N + -1)) = A^V-l)(/i + / 3 ) 2 

(JV_(JV_-1)) = A^V-l)(/ 2 + / 3 ) 2 



(N+N } = N(N-l)(fx + f 3 )(h + f 3 ) 

+Nf 3 . (80) 

The variance Vdyn, in the presence of resonances, is thus 
simply 



Vdyn — 



N(pi +P3)(P2 +P3) ' 



(81) 



One finds that the variance Vd yn increases with the 
fraction of resonances, p 3 produced in the final state. 
One also finds it to scale inversely to the number of par- 
ticles produced in the initial state. Note that in the limit 
p 3 = 0, Vdyn vanishes by our assumption of indepen- 
dent production. The simple treatment done here does 
not account for finite acceptance effects on the decay of 
resonances. Obviously, if too small a rapidity region is 
integrated, one of the decay partners may on average be 
missed, and \vd yn \ shall be increased accordingly. 

In AA collisions, one does not expect resonance pro- 
duction to be the sole cause of correlation, i.e. Vd yn < 0, 
but it is yet to be determined what fraction of the ob- 
served fluctuations may be attributed to resonance pro- 
duction or to truly dynamic correlations. In that respect, 
it shall be interesting to consider fluctuations of specific 
particle species such p/p in contrast to tt^ 1 or K ± given 
no known resonance decay into p + p whereas many res- 
onances exist that decay into 7r + + 7r~ or K + + K~ . 



X. 



SUMMARY AND CONCLUSIONS 



We introduced the net charge fluctuation measure, 
Vdyn, on the basis of two-particle correlation functions. 
We showed that for heavy ion collisions involving in- 
dependent nucleon collision and negligible rescattering 
of secondaries, Vdyn scales as the multiplicative inverse 
of the produced charged particle multiplicity. We also 
showed that Vd yn is simply related to other observables 
used or proposed for fluctuation measurement by various 
authors. We found however that the different fluctuation 
measures have different dependence on the experimen- 
tal acceptance, detector efficiency, and collision central- 
ity. We showed that v& yn has a weak dependence on the 
rapidity range used experimentally to measure the fluc- 
tuations provided the rapidity range is of the order or 
smaller than the two-particle correlator width, whereas 
observables such as $ have basically a linear dependence 
on the size of the acceptance. We found also that Vdyn 
is, by construction, independent, to first order, of the de- 
tection efficiency whereas measures such as $, luq have a 
explicit dependence on the detection efficiency. We also 
found that charge conservation has a finite, and actually 
sizable effect on the charge fluctuation measure Vdyn de- 
termined by the total charge particle multiplicity (over 
4-7T and independent of the detector acceptance used to 
measure the net charge fluctuations. We further showed, 
as also pointed out by Mrowczynski |l0| , that the $ 
measure shall be independent of the collision centrality 
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provided the collision dynamic is also independent of the 
collision centrality. Note however that because the detec- 
tion efficiency may be a subtle function of the detector 
occupancy, and hence the collision centrality, caution has 
to be exercised when interpreting uncorrected measure- 
ment of $ vs. collision centrality. Finally, we presented, 
as an example, a simple particle production model that 
can be used to account for the production of resonances 
as well as charge conservation. 
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The factor e a is the probability that a particle of type 
"a" is detected. One gets similarly for the 2nd moment 
and cross term: 

(»a) = p^Y, P D(™\M)P(M)[e 2 MM 2 + a 2 a M) 

^ ' M 

+£ a (l - e a )fi a M] , 
{<) = [e 2 a {°l - 1) + e«] (M) rn + e 2 y a (M 2 ) m , 



(A4) 
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(n a n b ) = e a e b fi a fi b ((M 2 ) m - {M) 2 m 

rn ■ 



(A5) 



APPENDIX A: FINITE EFFICIENCY EFFECTS 
ON THE MEASUREMENT OF v dyn 

A range of collision impact parameters is selected in 
experiments using a measured multiplicity m (or a sim- 
ilar observable). This introduces additional fluctuations 
because a single m corresponds to a range of impact pa- 
rameters. In this appendix we estimate the effect of cen- 
tralit y selection. We use these results in sees. [II and 

vtT| 



We assume, in the Gaussian approximation, that the 
moments scale with the true multiplicity M as 



(N a ) = VaM, 

(N 2 ) = £M 2 + olM, 
(N a N b ) = fia^M 2 + £ ab M 



(Al) 



where /i a and fj, b are average branching fractions for the 
production of species "a" and "b" respectively, while a 2 
and £ a t, are their variance and covariance. These relations 
are strictly true in the independent collision model or 
the wounded nucleon model where both M and iV a are 
respectively proportional to the number of sub-collisions 
or the number or strings. The first moment (|7l] ) is then 

(n a ) = - 5 ^Y / P D {m\M)P(M)e a n a M 
P(m) ^ 



?a»aj^Yl MP D (m\M)P{M) 



M 



e a fi a (M), 



(A2) 



where we have introduced the expectation value of M at 
fixed m defined as 

{M)m = Plm) MP D (m\M)P(M) (A3) 



The correlators R aa and R ab are therefore 

^a-Ma 1 , (M 2 ) m -(M) 2 m 



(MYt 



and 



Sab 1 (M 2 ) m -(M) 

A (M) m 



(M)'i 



(A6) 



(A7) 



The variance Vdyn — Raa + -^66 ~ 2Rab measured at a 
given m is then 



Vdyn{m) 



where 



v 



Ma 



(M) r 



M 



(A8) 



(A9) 



This expression amounts to the value of Vdyn evaluated 
at M = (M) m . One finds that the correlators R ab ex- 
hibit a contribution from the variance (M 2 ) m — (M)^ 
whose magnitude depends on the detector response func- 
tion width. The variance Vdyn however does not have 
such a contribution and as such is also independent of 
the detection efficiency for measuring M. 

Note that the above result implies that Vdyn is robust, 
i.e. independent of det ection efficiencies, in the Gaus- 
sian approximation (Al). An explicit dependence on ef- 
ficiencies would arise if the Gaussian approximation is 
not valid, e.g. if the detector response functions differ 
markedly from Binomial or Gaussian functions, or if the 
efficiencies exhibit very large variations with detector oc- 
cupency. 
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